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Abstract - We examine the force between two charged surfaces immersed in aqueous mixtures 
having a coexistence curve. For a homogeneous water-poor phase, as the distance between the 
surfaces is decreased, a water-rich phase condenses at a distance Dt in the range 1 — lOOnm. At 
this distance the osmotic pressure can become negative leading to a long-range attraction between 
the surfaces. The osmotic pressure vanishes at a distance < Dt, representing a very deep 
metastable or globally stable energetic state. We give analytical and numerical results for Dt and 
De on the Poisson-Boltzmann level. 



The forces between charged objects in electrolyte so- 
Ititions are of fundamental importance in biology and 
colloidal science. Within the Poisson-Boltzmann (PB) 
mean-field theory, in a single pure solvent, the inter- 
action between symmetrically charged colloids is always 
repulsive [l]. Experiments, on the other hand, have 
shown that highly charged colloids can attract each other 
when multivalent ions are present. This discrepancy has 
been explained theoretically by counterion correlations p] , 
fluctuation-induced forces [s] and other non-electrostatic 
interactions 'J. More recently, the PB theory was gener- 
alized by strong coupling theory (5][6], which predicts an 
attraction for large values of the so-called coupling param- 
eter [7l|8l. 



an important role not studied before 18 ■ 22 



In this Letter we use the PB theory to show that strong 
attractive forces appear between similarly charged colloids 
in mixtures. In mixtures one must also take into account 
that the medium itself becomes inhomogeneous due to di- 
electrophoretic and solvation-related forces. The preferen- 
tial solvation energy of ions in one of the solvents [l||9] is 
appreciable and can even be much larger than the thermal 
For a recent review on ion-specific solva- 
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tion effects within the PB theory see 
on preferential solvation in binary mixtures looked at the 
phase behavior in the bulk 14,15 , surface tension 16 



and the interaction between surfaces but not in immisci 
ble liquids 



17 



As is shown below, in partially miscible 
mixtures, the behavior is qualitatively different and strong 
forces occur. In colloidal suspensions these forces can have 



We consider two positively charged colloids in a mixture 
of polar solvents. The colloids are modeled as flat surfaces 
located a,t z = zLD/2 and uniformly charged with a charge 
density ea per unit area, where e is the elementary charge. 
A small amount of monovalent ions and weakly charged 
surfaces are assumed and therefore the coupling parameter 
is small and the PB theory is applicable [6] . The partially 
miscible solvents have a coexistence curve below the mix- 
ture critical temperature Tc in the absence of ions. This 
coexistence curve is further modified by the presence of 
ions 
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The grand potential density is given by 14 23 



+ ?i+ (log(uon+) - l) + (log(uon^) - l) 

- (Am+7i+ + Au"n")</)- A+n+ - A^n" -^0 (1) 

Here the Boltzmann constant is set to unity, T is the ther- 
mal energy and C is a positive constant, (p is the volume 
fraction of water (0 < (p < 1); far away from the sur- 
faces the mixture is homogeneous and water-poor with 
composition (po < 1/2 and ion densities uq. We use a 
regular solution form for the free energy of a binary mix- 
ture vofb = 01og((/)) + (1 - 0) log(l -(/.)+ - </>) p, 
where x ^ ^/T is the Flory parameter and Vq = cris 
the molecular volume of both liquids. The third term in 
Eq. ([T]) is the electrostatic energy of the mixture, where 
ip is the electrostatic potential. The mixture dielectric 
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constant depends on the composition through a hnear re- 
lation: = ec + {sw—£c)4>: where £c and the pure 
cosolvent and water dielectric constants, respectively. The 
second line of Eq. ([I]) is the ideal gas entropy of point-like 
ions, where n"'" and n~ are the positive and negative ion 
density, respectively. The first term on the third line is 
the bilinear solubility interaction of the ions and the sol- 
vent: the parameters Aw+ and Au~ measure the affinity 
of the positive and negative ions toward the water environ- 
ment, respectively. We use the common case where both 
ions are hydrophilic and assume the symmetric interac- 
tion Au~^ = Au~ — Au > 0. Lastly, and ^ are the 
Lagrange multipliers (chemical potentials) of the positive 
and negative ions and water composition, respectively. 
The ion densities obey the Boltzmann distribution 



(2) 



In a salt reservoir we have A* — \og{v()no) — Au^q. Alter- 
natively, when the mixture contains only negative counte- 
rions, A~ is determined self-consistently by charge conser- 
vation: Jn^dz — 2a. 

The electrostatic potential is determined by the Poisson 
equation 



V • (e((/))VV') = -e{n+ - n~) 



(3) 



supplemented by the boundary condition on the surfaces 
— n • V?A(±_D/2) — e(j/e{(p), where n is the outward unit 
normal to the surface. Finally, the Euler-Lagrange equa- 
tion for (h reads 



dfb 



1 de 



(VV')^ - AM(n+ + n") - ^ (4) 



In order to isolate solvation and electrostatic effects, we 
assume zero short-range chemical or long-range van der 
Waals interactions with the surfaces, leading to the bound- 
ary condition — n • V(/)(z = zLD/2) = 0. 

The net force exerted on the surfaces by the liquid is 
given by the osmotic pressure 11 = P^z — Pb, where Pb 
is the bulk pressure: Pb/T = 4>odfb{4>o)/d(j) - fb{4>o) + 
2no (1 — Au(f)o). —Pzz is zz component of the Maxwell 
stress tensor [l7l[25]: 
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dz 



,dfb 



-fb 



+ (I - AM(/»)(n+ + rT) 
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, de 
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di) 
dz 
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Pzz is independent of z and can be calculated at the mid- 
plane (z = 0) where by symmetry d(j)/dz — dtp/dz — 0. 
Multiplying Eq. Q by and inserting into Eq. ([s]) we 
obtain: 



n Tn.n 



' Tujh 



Pb 
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where uJb = fb — and (pm and the composition 

and total ion density at the midplane, respectively. 
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Fig. 1: (a) Composition profiles 4'{z) in the strong screening 
limit for D < Dt (dotted line), just before and after Dt (dashed 
and solid lines) and at De (dash-dot line). Here no — O.IM, 
Alt — 2 and a — lnm~^. (b) The same in the ideal gas regime 
with no salt, Au — 11 and a = O.Olnm"^. Here and, unless 
stated otherwise, in all other figures we took for the mixture 
(l>o = 0.4, T/Tc = 0.9915 and C = x/a [24]. As an approx- 
imation of a water-l-butanol mixture we used Tc — 398K, 
vo = 3x 10"^V^ Ea = 17.8 and e„ = 80. 



In aqueous mixtures, water is drawn towards the walls 
by field gradients (dielectrophoretic) and solvent-induced 

These forces are interde- 



(electrophoretic) forces 17 23 
pendent as is evident from Eqs. ([s]) and Q . From the 
solution of the governing equations we find two possible 
scenarios, distinguished according to the composition pro- 
files at infinite separation. In the first scenario, when the 
surfaces are far apart, (j){z) near each surface is composed 
of three layers: a thin water-rich layer (0 > 1/2) of width 
t at the surfaces, a second layer with width of the bulk 
correlation length ^ adjacent to it, where the composition 
4>^ (f>o decays to the bulk value, and a third "bulk" layer 
where ~ 0o. 

As the surface separation decreases down to D = Dt, 
the general features of this profile do not change. In partic- 
ular, only the width of the (f> (f>Q region changes if screen- 

— 1/2 

ing is strong, that is, when Id <^ D where ^ is 
the Debye length calculated at e — e{<j>o). The dashed 
curve in Fig. [ija) shows 0(z) just before Dt for this case. 

In the second scenario, dielectrophoretic and solvation 
forces are not strong enough and at large surface separa- 
tions <j) is always close to the bulk value (j)o. As D de- 
creases from infinity the water composition between the 
surfaces increases continuously. In particular, in the ideal 
gas regime, defined hy b ^ D, where b oc cr~^ is the Gouy- 
Chapman length, 4>{z) is nearly uniform. Fig. [ijb) shows 
such profiles for a mixture containing only counterions. 
Here, the ion profile n{z) and electric field are nearly uni- 
form leading to small composition gradients. 

In both scenarios, when D < Dt the whole space be- 
tween the surfaces becomes rich in water (see the solid 
curves in Fig. [T]). The transition occurs at a distance 
in the range Dt = 1-lOOnm depending on the parameter 
values, and is accompanied by a decrease in the osmotic 
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Fig. 2: (a) Osmotic pressure 11 vs surface separation D in 
tiie strong screening limit for different values of no [other 
parameters as in Fig. [ija)]. Inset: excess grand poten- 
tial for no = O.IM. solid (dashed) line corresponds to stable 
(metastable) solutions, (b) IT in the ideal gas regime for differ- 
ent values of Au [other parameters as in Fig. [ijb)]. Inset; the 
function —ujb{(j>). Filled square and circle markers correspond 
to 00 and the binodal compositions, respectively. In (a) and 
(b) open markers on solid lines correspond to distances D with 
the same symbols in Fig. [l] (a) and (b) . 

pressure. This can lead to a negative osmotic pressure, 
such that a surface separation D,, < Dt exists at which 
n = 0. (j){z) at D = De is shown in the dash-dot curves 
in Fig. [1} 

The excess grand potential per unit area relative to infi- 
nite separation is An{D) = - J^n{D')dD' . An{D) has 
a cusp at D = Dt and this results in a discontinuity of 
the pressure. Close to Dt, both a water-poor and a water- 
rich profiles are solutions of the Euler-Lagrange equations 
but only one of these is a stable solution. At D > Dt 
the water-poor phase is stable, at D < Dt the water-rich 
phase is stable, and at D = Dt the grand potentials of the 
two phases are equal. Furthermore, Afl has a minimum at 
Df, (De < Dt) corresponding to mechanical equilibrium. 
An example of AVl(D) is shown in the inset of Fig. [2ja), 
where the dashed line is AH. of metastable solutions. 

The transition to a water-rich phase has a different phys- 




Fig. 3: Osmotic pressure in the weak screening regime with 
no = 10~''M and a = lnm~^ for three values of Au. Inset: the 
corresponding excess grand potential AQ. 

ical origin in the two limiting cases. In the strong screening 
regime, the transition is promoted by the energy gained 
when the interface vanishes, as in capillary condensation. 
The thin water-rich layer near the walls remains nearly 
independent on D. In the ideal gas regime, on the other 
hand, when the composition approaches the coexistence 
composition, preferential solvation promotes a transition 
to a water-rich phase. For intermediate cases, both mech- 
anisms play a role in the transition. These effects are 
enhanced but are not limited to the vicinity of the criti- 
cal temperature where the differences between the phases 
become smaller. 

In Fig. [2] we plot the osmotic pressure as a function 
of surface spacing; when 11 > the surfaces repel each 
other while for 11 < they attract. 11 is discontinuous at 
D^Dt. InFig.[2l;a) we show n for different values of no 
in the strong screening limit. An increase in no decreases 
n at large distances but increases it at small distances 
(entropy loss of the ions). 

The negative jump in II at D — Dt increases with de- 
creasing no- This can be explained by the interplay be- 
tween the first two terms in Eq. At D — Dt the in- 
crease in (f)m dominates the decrease in 11, an effect more 
pronounced for smaller values of no since the positive ideal 
gas term is proportional to n^ cx uq. Recall that in a pure 
solvent only the ideal gas term exists and hence 11 is al- 
ways positive. Furthermore, this entropic repulsion will 
eventually dominate 11 leading to 11 = at a separation 

De. 

In Fig. [2jb) we plot 11 for different values of Am in 
the ideal gas regime. For Am = 1 the pressure is purely 
repulsive since preferential solvation is not strong enough 
to induce a water-rich phase. For Au = 7 the interaction 
is repulsive down to Dt where 11 becomes negative. When 
Au increases to 11 the transition is at a larger distance 
and to a more negative pressure. The inset of Fig. [2jb) 



p-3 



S. Samin et al. 



shows the function — ajf,(0), showing a decrease in pressure 
at the transition [cf. Eq. 

Fig. [Sjshows n in the weak screening regime where Id ~ 
D and for three different values of Aw. When Am = 0, the 
dielectrophoretic force alone can induce a water-rich phase 
albeit the pressure is always repulsive. The pressure can 
become attractive for Au = 1 or 7. Corresponding curves 
of Afl are plotted in the inset of Fig. [3j These show a 
metastable minimum at a finite value oi D, D = D^,, for 
Am = 1 and a global minimum for Am = 7. The depth of 
the minimum for Am = 7 is « 440T for two colloids with 
an effective surface area of lOOnmxlOOnm. 

For Au = 10, a second metastable minimum can ap- 
pear in the curve A^l{D) at a smaller value of D, see the 
dashed and dash-dot curves in Fig. |4j A similar value of 
Am is cited in the literature for mixtures of water and 1- 
butanol containing NaCl at room temperature 12 . When 
(7 is further increased (solid curves), this minimum can be- 
come globally stable. For large enough a, only the second 
minimum exists (not shown). In this large Ait case, pref- 
erential solvation leads to the liquid between the surfaces 
being nearly pure water; 0^ is close to unity and Wf, di- 
verges. 

The distance Dt can be obtained in the two limiting 
cases presented in Fig. [T] In the strong screening limit we 
use continuity of A51 at Dt to obtain [26|: 



D/2- 



aj((/>(z))dz 



Ph - Pi 



h 



(7) 



Ph, the pressure in the water-rich phase, 
from Ph = T{n+ + n^) - Tujb{^h) + Pb, 

from 



is calculated 



where 



+ 



are the midplane values and (jih is obtained 



dfb 



{(t>h) - 2AMnoe'^"('^''^'^°) - Mo = 



(8) 



where we used sa at the midplane. 

At D = Dt the surfaces can be regarded as isolated 
and a zeroth-order approximation of can be obtained 
using the well known result for a single surface with 
a homogeneous (p k (ph and a modified Debye length 
(e(<^ft)T/(2e2noe^"(^'--^«)))'^'. Since at \z\ < D/2-t 
eip/T <^ Au{<j)h — 0o), t is obtained from the condition 
eip/T = 0.01Am((/)/j — <Pq). In analogy with classical mean 
field theory [24] , the Landau expansion of Eq. ([s]) around 
(pc = 1/2 gives for ^ 
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(9) 



The preferential solvation term in the root is comparable 
in magnitude to 1 — T/T^. Thus, the bulk correlation 
length is modified appreciably by the preferential solvation 
of the ions. 

Fig. [Sf^a) shows the comparison of Dt from Eq. Q with 
numerical results for different values of no. The agreement 
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Fig. 4: (a) osmotic pressure for a large Au — 10, no = 10~*M 
and three values of a. (b) the corresponding excess grand po- 
tential. For a = 0.54nm~^, AO, has only one minimum (dash- 
dot curve) at D < Dt. A second, metastable minimum appears 
for a — 0.56nm~^ (dashed curve), whereas for a — 0.62nm~^ 
(solid curve) this minimum becomes globally stable. 
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Fig. 5: Transition and equilibrium separations Dt (circles) and 
De (squares), (a) The strong screening case vs no. (b) The 
ideal gas regime vs a. Solid curves are analytical expressions 
given in the text. In (a) we approximated (j>{z) by a linear 
decrease from ph to po at the interface region. In (b), for a 
beyond the dashed curve, a water-rich layer near the surfaces 
exists even as D — >■ oo. 



is quite good despite the crude approximation. D,, in Fig. 

decreases slightly when uq increases, its value being 
w 5nm. 

In the ideal gas regime p{z) ~ const, and the counter 
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Fig. 6: Interaction potential for a mixture with a critical com- 
position (f)o = 1/2 at three temperatures above Tc- For t = 1.01 
(dash-dot curve) the interaction is purely repulsive. When 
the temperature is decreased the interaction turns attractive 
(dashed curve), becoming stronger closer to Tc (solid curve). 
Here no = O.OIM, Am = 4 and a — Inm^^. The walls have a 
surface oi A — 0.01/im^. 



ions density 
density 
reads 



24 



' is uniform and equal to the average charge 
n^) ~ 2a/D. The composition equation 



dft 2Aua 



- A*o = . 



(10) 



Here preferential solvation merely shifts the chemical po- 
tential of the mixture. Hence, Dt occurs when is the 
binodal composition at which dfb/d4> = 0. Thus, 



(11) 



which means that Z?f — > cxd when the bulk composition 
approaches the binodal (/iq — > 0). 

De can be found by noting that H in Eq. ([6| gives 
Tie — uJbi(l>e) + Pb/T, whcrc (j)e and Ue are the composition 
and ion density at D^, respectively. Inserting this into Eq. 
^ we obtain for (j),. 



dtoh 



Au{ub{c^e)+PblT) 



(12) 



which is solved and then used to get rig and D^, from the 
equations for H and (n~), respectively. Comparison of 
the formulae for Dt and D^. with the full numerical results 
are presented in Fig. [5jb). As expected, the agreement 
is good for small values of a where the ideal gas limit is 
valid. The dashed curve in Fig. [sjb) marks the charge 
density above which a water-rich layer exists at infinite 
separation and the above approximations no longer hold. 

We find that at an attractive interaction is also pos- 
sible at T > Tc- In this case, since there is no misci- 
bility gap above Tc, the water-reach transition is missing 



and hence no discontinuity in the pressure is observed. 
Instead the pressure becomes smoothly attractive as the 
plate distance is reduced. Fig. |6] shows the interaction 
potential U — Af2 x A for two plates in a critical mixture 
at three different temperatures. As the critical tempera- 
ture is approached, a purely repulsive potential (dash-dot 
curve) becomes attractive (dashed curve). The interaction 
is more attractive closer to Tc (solid curve). The interac- 
tion energy and length scales and the temperature trend 
shown in Fig. [6] are similar to those observed in recent ex- 
periments on the salt-dependent interaction of a charged 
particle suspended in a critical binary liquid mixture near 
a charged wall ^21^. In light of this, we believe that the 
mechanism we describe may be important to capture cor- 
rectly electrostatic effects in these experiments. Thus, it 
is possible to attribute some of the effects shown in these 



experiments to solvation related forces 27 in addition to 
the critical Casimir forces. 

The qualitative difference in the interaction between 
charged colloids or macromolecules in mixtures compared 
to pure solvents stems from the contribution to H of the 
second term in Eq. (|6]): u)h{(j)ra) ~ the midplane mixture 
grand potential. This term is absent in pure solvents. 
Thus, unlike similar ion induced phase transitions in pure 
solvents [28 we observe a jump to a negative pressure (at 
Dt ~ 1 — lOOnm). Moreover, the nontrivial dependence 
of ) on the system parameters {T,(j)o,no,a) through 

the governing equations leads to qualitatively different be- 
havior compared to usual condensation transitions due to 
surface fields 
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In addition the attractive force is in 
many cases strong and long-range compared to the van 
der Waals forc e [l] as is evident by the values of H at 
D = Dt [Fig. [2| and Fig. [S] and the energy minimum 
being deeper than ~ lOOT. 

The mechanism we describe should be at play in the 
aggregation of charged colloids in mixtures near the co- 
existence temperature 18 19 . We believe it is directly 



relevant to the attraction seen between colloids and sur- 
faces in mixtures also far from Tc and attributed to critical 
Casimir forces 
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